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ABSTRACT 

Our multi-epoch survey of ~ 20 sq. deg. of the Canis Major overdensity has 
detected only 10 RR Lyrae stars (RRLS). We show that this number is consistent 
with the number expected from the Galactic halo and thick disk populations 
alone, leaving no excess that can be attributed to the dwarf spheroidal (dSph) 
galaxy that some authors have proposed as the origin of the CMa overdensity. 
If this galaxy resembles the dSph satellites of the Milky Way and of M31 and 
has the putative My ~ —14.5, our survey should have detected several tens of 
RRLS. Even if My < —12, the expected excess is > 10, which is not observed. 
Either the old stellar population of this galaxy has unique properties or, as others 
have argued before, the CMa overdensity is produced by the thin and thick disk 
and spiral arm populations of the Milky Way and not by a collision with a dSph 
satellite galaxy. 
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Subject headings: stars: variables: other, Galaxy: stellar content, Galaxy: struc- 
ture 



INTRODUCTION 



The Canis Major overdensity was discovered by iMartin et al.l (j2004al ) as an excess in 
the density of 2MASS M-giant stars in the direction (/, b) = (240°, —8°) compared to similar 
fields in the northern galactic hemisphere. It is estimated to span ~ 100 sq.deg. on the 
sky and to be concen trated at a distance of ~ 7.2 kpc from the Sun (IMartin et al.l l2004bl : 
Bellazzini et al.ll2006l ). The large numbers of M-giant and red clump (RC) s tars in the CMa 
overd ensity suggest that its main population is of intermediate age ~ 7 Gyr (IBellazzini et al. 



2004T). A "Blue Plume" of stars in the color-magnitude diagram (CMD) has been interpreted 
by IBellazzini et al.l (120041 ) and iMartmez-Delgado et al.l (120051 ) as the main-sequence of a 
younger population (age < lGyr) 



The CMa overdensity was interpreted by IMartin et al.l (j2004al ) as the remnant core of 
a disrupting dwarf galaxy. However, its origin remains a topic of considerable debate. The 
dwarf galaxy hypothesis has been contested by several authors who, using different tracers 
and techniques, have proposed t hat the overdensity is produced by t he stellar population s 
of the warped Galactic thin disk (iMomanv et al.ll2004T). the thick disk dCarraro et al.ll2008l ). 



and the outer Norma-Cygnus sp iral arm dMoitinho et al. 



2006: Carraro et al. 



20051) . Cou nter 



arguments have been made by IMartin et al.l (j2004bl . 120051 ) and IBellazzini et al.l (120061 ) , in 
support of the dwarf galaxy nature of CMa, claiming that the Galactic warp is not sufficient 
to account for the overdensity. The pro per motions of a sample of blue plume stars have 
been measured by iDinescu et al.l (120051 ) . who concluded that these stars do not have the 
space motions e xpected of galactic sta rs, if they have the same radial velocity as the sample 
of M-giants that IMartin et al.l (j2004bl ) attributed to the CMa dwarf galaxy. 



The lack of consensus on this issue is probably a consequence of our limited knowledge 
of the complex Galactic structure in this direction and also the large and highly variable 
interstellar extinction. CMa lies in the third Galactic quadrant where the warp of the thin 
disk has its maximum amplitude (at I ~ 270°, IMomanv et "aDl2004l: lL6pez-Corredoirall2006l ) 
and where the Galactic thick disk might be warped as well (jCarraro et al.ll2008l). Also, in this 



quadrant the spiral arms have not been traced as well as in other directions (jCarraro et al. 
20051 : iMcClure-Griffiths et al.ll2005h . 



In Table HI we summarize the previous studies of the stellar populations in the CMa 
overdensity. Most are photometric and use techniques such as main sequence (MS) fitting 
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and modelling of either Hess diagrams (HD) or CMDs to obtain the age, me t allicity and dis- 



tance for the stellar populat ion of t he CMa overdensity . ICarraro et al.l (120051 ) ; lMoitinho et al. 



tood ): ICarraro et al.l (boost ) and IPowell etaD (booa i have made use of multiband photom- 
etry and two color diagrams (TCD) to reduce the effects of reddening. The top part of the 
table shows the studies of the young population in the region, while the bottom part indi- 
cates the intermediate age ones. As one can see from the table, no old population (> lOGyr) 
has been identified. 

In summary, there are two conflicting views of the nature of the CMa overdensity: 



The CMa overdensity is a dSph galaxy undergoing tidal disruption. The stellar pop- 
ulation of this galaxy is dominated by an int ermediate-age popula ti on, and there is 
a minority population of much you nger stars (IMartin et all l2004al B 1 Bel lazzini et al.l 
20061 : iMartmez-Dekado et allbood ). 



The overdensity in the inter mediate-age populat i on is produced by a projection of the 
warp ed thin and thick disks (IMomany et all2004lL6pez-Corredoirall2006l : ICarraro et al. 
20081 ) into the line of sight. The young population corresponds to the outer part of 
the Norma-Cygnus spiral arm which, being part of the thin disk, also follows the warp 
JCarraro et al.lbood : iMoitinho et al.l l20od : IConn et al.l 120071 ) . 



The goal of our work was to conduct a large scale survey of RR Lyrae stars (RRLS) in 
order to trace the oldest stellar population (> lOGyr)0. Old stellar populations containing 
RRLS have been identified in every d Sph satellite of the Milky Way (MW) and of M31 that 
has been properly surveyed (see e.g. I Vivas fc Zinnll2006h . If the CMa galaxy is similar to 
these systems in its stellar mix and has the proposed luminosity, a large excess (> 100) of 
RRLS is expected (see §4). Because RRLS are excellent standard candles, the distances to 
individual RRLS can be obt ained with high prec ision (< 10%). If there is a dSph galaxy 
at a distance of ~ 7.2 kpc (IMartin et al.l l2004al ) . an overdensity of RRLS is expected at 
magnitudes V ~ 15. The detection of an overdensity of RRLS that spans a narrow range of 
magnitudes would be conclusive proof that a dSph galaxy produces the CMa overdensity. On 
the other hand, the absence of an excess of RRLS would weaken the dwarf galaxy hypothesis, 
leaving more likely a Galactic origin for the CMa overdensity. 



Earlier progress reports of our survey can be found in iMateu et al.l (fcood . l2007h 
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PHOTOMETRIC OBSERVATIONS 



The large field of view (2?2 x 2?3) provided by the QUEST-I camera (IBaltav et al.ll2002h . 
installed on the lm Jurgen Stock Telescope at the Venezuelan National Observatory of Llano 
del Hato, made this an ideal instrument to conduct a large-scale, multi-epoch survey of an 
extended structure such as the CMa overdensity. The observations were made in the V band 
and were obtained in two runs during 2004 and 2005, centered at I ±2 240°, b ^ —9.5° and 
spanning ~ 30 sq. degrees with an effective coverage of ~ 20sq. deg. (which takes into 
account the area missed because of the gaps between the CCDs in the mosaic camera and 
the saturation of very bright stars). On average, we obtained ~ 19 epochs per field in 24 
nights, with a time spacing between consecutive observations ranging from a few hours up 
to several weeks. Infrared J, H and K magnitudes from the 2MASS catalogue were also 
obtained for each star. 

The area covered by our survey and previous studies is shown in Fig. EH where a 
schematic diagram is overlaid on an extinction map derived from the ISchlegel et al.l (119981 ) 
dust maps. The CMa feature spans ~ 100 sq. deg. and has a smoo thly varying density 
across the plane of the sky, which according to iMartin et al.l (|2004al ) can be fitted with 
Gaussian profiles. The axes of the ellipse shown in Fig{T] correspond to the full-width at 
half maxima (FWHM) derived for CMa using th e aspe ct ratio of ~ 5 : 1 and the FWHM in 
the latitude direction measured by iButler et al.l (120071 ). Assuming a Gaussian profile with 
these FWHM we computed that ~ 3% of the stars in CMa overdensity shoul d be present 



i n our survey's footprint. The percentage could be smaller than this because IButler et al. 



(2007) considered 5 : 1 to be a lower limit of the aspect ratio. From Fig. 6 of IMartin et al. 
( I2004al ). we estimate that the aspect ratio is probably not much larger than ~ 6 : 1, and 



consequently, that our survey's coverage is > 2%. In the following discussion, we adopt 2% 
to be conservative. As shown in Fig. [IJ our survey includes the same central regions where 
others (e.g., Martinez-Delgado et al. 2005; Bellazzini et al.2006; Butler et al. 2007) have 
detected large overdensities of stars. One can also see that the extinction is highly variable 
in the vi cinity of CMa and that i ts main body is located in a low extinction window, as 
noted by iRocha-Pinto et al. (bood ). 
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2.1. Image Reduction, Photometry and Astrometry 

The image frames were reduced using standard iTMi^ routines for bias subtraction and 
flat-fielding. Due to the low galactic latitude of the surveyed region, the fields are highly 
crowded, thus requiring the computation of magnitudes through PSF photometry. The PSF 
photometry was obtained using DAOPHOT tasks with a model PSF spatially variable to 
second order, constructed from ~ 130 — 150 stars per image. 



The astrometric solution was computed for each image using the program CM1 ([Stock 



19811 ). which calculates the astrometric matrix based on t he coordinates of y 150 catalogue 



stars per image, identified from the USNOB1.0 catalogue ( Monet et al.ll2003l ). The precision 
of the obtained astrometric solutions was 0'.'18, which was calculated by matching the coor- 
dinates of all objects with the 2MASS point source catalogue, and corresponds to the mean 
difference between the CM1 and the 2MASS coordinates for each object. This precision is 
sufficient for identifying each star in the different observations. 



2.2. Photometric Calibration and Extinction Correction 

The photometric calibration was made in two steps, the normalization of the instrumen- 
tal magnitudes with respect to a reference night, and the calibration of the reference night 
photometry using secondary photometric standards. 

In the first step, the magnitude zero-point of each night with respect to the reference 
night was obtained by calculating a 3-sigma clipped mean of the magnitude differences of 
the stars that were detected on both nights. The use of ~ 2000 stars in the calculations 
made these zero-points statistically robust. The obtained zero-point was then added to the 
instrumental magnitudes of all the objects observed on that particular night. At the comple- 
tion of this process, mean photometric errors were evaluated as a function of magnitude by 
calculating the standard deviation of the multi-epoch observations of each object. A typical 
plot of standard deviation a of PSF magnitudes as a function of the mean magnitude is 
shown in Fig. [2l 

For the second step, observations of the secondary standard fields were made with the 
1.3m telescope of the SMARTS consortium at the Cerro Tololo Interamerican Observatory 
(CTIO), Chile, with the 2048 x 2048 pixel Fairchild 447 CCD camera, which has a field 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asoci- 
ation of Universities for Research in Astronomy Inc., under cooperative agreement with the National Science 
Foundation 
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of view of 6' x 6'. The observations were made in the R and V filters on 2004 January 
23, under photometric conditions. A field containing Landolt standards was also observed 
during the night at different airmasses. The photometric solution for the night was applied 
to the objects detected in the secondary standard fields. The 316 stars that were selected 
as secondary standards have photometry in both filters and photometric errors < 0.015 
magnitudes. 

The secondary standards fields were chosen to be centered on open clusters lying in 
the region, to ensure having blue stars among the secondary standards. Consequently, the 
color distribution of the 316 secondary standards was quite broad, 0.1 < V — R < 0.8. Each 
row of CCDs in the QUEST camera was calibrated independently. The mean rms of the 
calculated solutions is 0.03 magnitudes. 

The limiting magnitude obtained is Vu m = 19.5, defined as the typical magnitude of 
objects with o = 0.1 mag. The completeness magnitude obtained is V com = 18.5, which 
corresponds to the mode of the magnitude distribution of all objects. Finally, the saturation 
magnitude of the survey is V sat ~ 13.5. The heliocentric distance range explored by our 
RRLS survey is therefore 3 kpc < D & < 39 kpc (for the complete sample and assuming 
My = +0.55), which clearly includes the distance range of the CMa overdensity (5 kpc 
< D & < 10 kpc, see Tabled]). 

Finally, the extinction correction was performed on a star by star basis, since extinction 
is highly variable in the surveyed region (as shown in Fig. []JJ due to its proximity to the 
Galactic plane. This was done by subtracting the corresponding Ay at the position of 



each star, computed from a linear interpolation of the ISchlegel et al.l (119981) dust maps, and 
including the asymptotic correction proposed by Bonifacio et ali ( 200ol ). 



3. RR LYRAE STARS 



3.1. RR Lyrae search 



The first step in the RRLS search was the identification of s tars showing photo metric 
variability in the V magnitude, using the procedure explained in I Vivas et al.l (12004 ). This 
consisted of computing, via a \ 2 test, the probability of the V magnitude dispersion being 
solely due to the photometric errors. From the catalogue of 286, 837 stars, 2775 stars had 
P x 2 < 0.001 which indicates that they have a probability of 99.9% of being variable. 



The RRLS search was further restricted to the 651 variable candidates fulfilling the 
following criteria: (i) intrinsic color indices (V — J)q, (V — H) , (V — K) consistent with 
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those of spectral types earlier than GO. (ii) observed amplitude AV > 0.20, and (iii) number 
of observations in the V filter, n\ 



The RRLS search followed the procedure devised by iLaydenl (119981 ) for fitting the light 
curves of RRLS. This consisted of fitting light-curve templates while varying four parameters: 
period, amplitude, phase offset and magnit ude at maximu m light. For the fits of RRab 
variables, we used the six templates made by ILaydenl ( 119981 ) . which cover the range of light- 
curve shapes of this type of star. For RRc stars, we used either a cosine curve or a template 
made by averaging and smoo thing the light cu rves of 10 well-observed RRc stars from the 
QUEST RR Lyrae catalogue J Vivas et al.lbooj ). 



The ranges for the trial parameters were expl ored according to the known distribution 
of periods and amplitudes for RRLS of each type (jVivas et al.ll2004l ). using a step of 10 -5 d 
for the period in the range 0.4 — 1.2 d for RRafr and 0.2 — 0.5 d for RRc; a step of 0.01 mag 
for the amplitude in the ranges 0.2 — 1.4 mag for RRab and 0.18 — 0.7 for RRc; and steps of 
0.01 mag and 0.005 for the V magnitude at maximum light and the phase offset respectively. 
These values were then taken as the uncertainties of the corresponding parameters. For 
each set of trial parameters, the \ 2 °f the fit was computed and finally, the four best fits 
for each RRLS type were visually inspected for every star. The stars for which an adequate 
template fit was found were selected as initial RRLS candidates. Additionally, the epoch 
at maximum light and the flux-averaged V magnitude (V) were calculated using the best 
fitting light-curve parameters. 



3.2. Photometric Follow-up 

The observations with the Stock telescope yielded for some RRLS candidates poorly 
sampled light curves. Consequently, we obtained additional observations with the 0.9m, 1.0m 
and 1.3m telescopes of the SMARTS consortium at Cerro Tololo Interamerican Observatory 
(CTIO), during Sept-Oct 2004 and Jan- Feb 2008; and with the 1.0m Reflector at the National 
Observatory of Venezuela (NOV), during Jan 2008. All observations were made in the V 
filter. This photometry was normalized to the reference night of the QUEST photometry 
using ~ 120 — 250 field stars per image for the CTIO observations and ~ 70 — 80 for the 
NOV observations. On average, 5 additional epochs were obtained for each candidate RRLS, 
yielding an average 24 observations per star. Light curve parameters were recalculated with 
the additional observations, including exploration of period aliases. With these additional 
observations, only 16 stars passed our tests for RRLS candidates, 10 of type c and 6 of type 
ab. 
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3.3. Spectroscopic Follow-up 



Spectroscopic follow-up of the 14 brighther (V < 16.5, equivalent to D & < 12 kpc) 
RRLS was conducted in order to obtain radial velocities and metallicities. These bright 
stars lie within the distance range of the CMa overdensity. 

Spectra were obtained with the Cassegrein spectrograph at the 1.5m SMARTS telescope 
in Cerro Tololo Interamerican Observatory, Chile. The observations were obtained in ser- 
vice mode in the months 2004 November-December, 2005 January, and 2008 February. The 
instrumental setup was the same used by our gr oup to obser v e a la rge number of halo RRLS 
from the QUEST survey. We refer the reader to IVivas et al.l (120081 ) for details on the instru- 
mental setup, data reduction, and methods to estimate radial velocities and metallicities. 
Here we will summarize only the more important issues. The spectra have a resolution of 
4.3 A and a spectral range from 3500 to 5300 A. Radial velocities were obtained by Fourier 
cross-correlation (IRAF's fxcor task) with several radial velocity stand ards of the sam e spec- 
tral type. The radial velocity standards were chosen from the list of iLaydenl (119941 ) which 
also serve as pseudo-equivalent width standards for estimating the metallicities. Each RRLS 
was observed twice at different phases during the pulsation cycle. To obtain the systemic 
velocity, V 1 , of each star, t he radial velocity curve of the well-studied star X Arietis was 
fitted to the data (see also IVivas et al.l 120051 ) . The metallicities were estimated using the 



calibration by ILaydenl (119941 ) of the pseudo-equivalent width of the Ca II K line (W(K)), 
corrected for interstellar absorption, and the mean pseudo-equivalent widths of the /3, 7, and 
S Balmer lines (W( H)). The estimated error of the [Fe/H] measurements, which are on the 



Zinn fc Westl ( 119841 ) metallicity scale, is 0.15 dex for the RRafe and 0.20 dex for the RRc. 



3.4. Contaminants of the Sample 



Our sample is expected to be contaminated by other types of variable stars. In order 
of decreasing likelihood, these are contact eclipsing binaries (W UMa stars ) of A-F spectral 
type, pulsating blue straggler stars of the Delta Scuti and SX Phoenicis types, and Cepheids 
of type II and of the anomalous classification. 



In many cases the light curves of W UMa stars resemble a cosine curve (lRucinskilll993l ). 
making it difficult to distinguish them photometrically from RRc stars, since their period 
ranges also overlap. However, several orbital configurations do have light curves that can be 
distinguished, and consequently, we also fitted a W UMa templates to some of the candidate 
RRc stars. This template was co nstructed from the Fourier coefficientes of iRucinskil (119931 ) 
(i = 80°, q = OA J = 0), following ILaydenl (I1998I ). We used only this W UMa template, which 
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we consider representative of the type, since our goal was only to reduce the contamination 
that these stars introduce in the RRc sample, and not to survey their properties. The W 
UMa template yielded a better fit than the RRc template for 3 stars, namely 25711, 10651 
and 18056. The radial velocities of these stars are also consistent with the interpretation 
that these stars are not RRLS. The radial velocities of stars 25711 and 10651 appear to be 
constant, while for star 18056 the radial velocity curve template for RRc stars does not fit 
the observations properly. Furthermore, the mean velocities of these 3 stars are 81 km s _1 , 
26 km s _1 and 77 km s _1 , respectively, which lie close to the radial velocity peak expected 
for thin disk stars, which for the survey region is Vy ~ 69 km s _1 with aV 7 ~ 22 km s" 1 
(see §4.2p . We conclude that stars 25711, 10651 and 18056 are most likely eclipsing binaries 
from the thin disk, and have classified them as W UMa stars. 

In the selection of RRLS candidates, we imposed color cuts to exclude the numerous 
red variables in the fields, but some non-RRLS may pass these criteria. Our spectroscopy 
revealed that the RRc candidates 6460, 5564, and 17237 have late F spectral types. Since 
these types are too late for RRc stars (Smith 1995), we have removed them from the sample 
of candidate RRLS. 

The Delta Scuti and SX Phoenicis variables, also known as dwarf cepheids, have spectral 
types A0-F5 III-V, V amplitudes < 0.7 mag and periods in the ranges 0.01 — 0.2 d and 
0.04 — 0.08 d respectively. Consequently, the period ran ge of Delta Scuti stars overlaps wit h 



the short period end of the RRc period distribution (iSmithl 119951 ; iKholopov et al.l 119981 ). 
In the case of SX Phoenicis stars, it is expected that some stars will have period aliases 
within the range expected for RRc stars, which may escape detection by the low number 
of epochs and the time spacing of our observations. Although much less numerous than 
eclipsing binaries, Delta Scuti and SX Phoenicis stars are expected in the Galactic disk and 
halo respectively, and hence are possible contaminants of the RRc sample. 

Finally, a possible contaminant of the KRab sample are Anomalous Cepheids and short 



period type II Cepheids, als o known as BL Herculis and W Virginis stars (ISandage &: Tammann 



20061 ; IKholopov et al.l ll998). These stars have light curve shapes and period ranges that par- 
tially overlap with RRLS stars. However, these stars are very rare in comparison to the RRLS 
in the halo population, which is the only MW population in which they have been found. 
Given that the halo population is itself a small contributor to the stellar density in the CMa 
fields (see §5), the contamination from these stars should be small compared to the other 
types of stars. 

After the removal of the likely contaminants, 10 stars remained in our sample of RRLS 
candidates. The light curves for the 10 RRLS are shown in Figure [31 and those of the 3 
W UMa stars are shown in Figure HI The light curve and physical parameters obtained 
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for the 10 RRLS in Table [2], including heliocentric distances of RR LS calculated based on 
the flux-averaged mean magnitude (Vo) and assuming My R = +0.55 flDemarque et al.ll2000l : 



Vivas et al.l l200lh . The uncer tainty in the distance to the RRLS has been estimated in 7%, 



following I Vivas & Zinnl (120061 ). The light curve parameters of the rejected RRLS candidates 
are summarized in Table [31 We must caution, however, that some contamination by non- 
RRLS may be still present, particularly for the sample of RRc stars. 



3.5. Completeness of the Survey 

The survey's completeness was evaluated as the fraction of synthetically generated light 
curves for which the period was succesfully recovered. A total of 3600 synthetic RRLS 
light curves were generated for both RRa6 and RRc types, with 10, 15 and 20 data points 
and a time sampling characteristic of the survey. Typical random photometric errors were 
added as a function of magnitude (Fig. [2]). The synthetic light curves were used as input 
for the RRLS identification algorithm described in §3.11 If for a synthetic light curve, the 
difference between the real and the best fitting period was better than 1%, it was considered 
as recovered. The percentage of recovered light curves as a function of magnitude for types 
RRafe and RRc are shown in Fig. [51 For this plot, we averaged the results obtained for light 
curves with different numbers of data points. One can see from this figure that the average 
completeness of the survey is ~ 90%. 

RRLS 26896 has been previo usly identified as an RRLS by the All Sky Automated 



Survey (ASAS-3. lPojmanskill2002l ). The period derived for this star in the ASAS-3 survey 
(Pasas — 0.39725 d) coincides with our best-fitting period (see Table [2]). The ASAS V- 
filter observations are systematically brighter than ours, but this is due to a bright neighbor 
(V ~ 13.6) of star RRLS 26896, lying at a distance of ~ 14", which is unresolved in the 
ASAS survey. This is the only RRLS from the GCVS or the ASAS-3 that lies in the surveyed 
region, and it was successfully recovered by our survey. 



4. EXPECTED NUMBER AND DISTRIBUTION OF GALACTIC RRLS 

RRLS are found in the halo and thick disk stellar populations of the MW, and a certain 
number of them are expected in the region that we have surveyed. Below we estimate the 
numbers of these stars and their distributions in space and radial velocity to see if there is 
an excess of RRLS that should be attributed to another source. 
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4.1. Expected Radial Profile of Galactic RRLS 

Since the volume of our survey is naturally expressed in heliocentric Galactic coordinates 
(D Q , I, b), the number Nrr of RRLS in a given volume is expressed as 



N, 



RR 



D f 



I)" 

Q 



b Jl 



p(D Q , I, &)-Dq dD Q cos b db dl 



This equation requires that the density profile p be expressed in terms of heliocentric 
coordinates. Since the density profiles of the halo and thick disk are naturally written in 
galactocentric coordinates, the coordinate transformation complicates the integral in Eq{TJ 
making extremely difficult to obtain an analytical solution. Consequently, we solved it nu- 
merically by using a Monte Carlo implementation of the Von Neumann Rejection Technique 



flPress et al. 1120021 ). 



The R RLS density profile of the Galactic halo has b een well constrained by several 
authors (e.g. lVivas fc Zinnll2006l ; IWetterer fc McGrawlll996l ). In the case of the Galactic thick 
disk, the RRLS densi ty profile has only been characterized in the direction perpendicular to 
the plane of the disk (jLaydenlll995l ). which is due to the difficulties in surveying the crowded 
fields near the Galactic plane where the interstellar extinction is high. The density profile 
of t he thick disk has been studied by several authors using different tracers a nd techniques 
(e.g. iBrown et al.ll2008l : ICabrera-Lavers et al.ll2005l ; lLarsen fc Humphreysll2003l ) , and we have 
adopted these profiles as valid for the RRLS too. The assumed analytical expressions for the 
density profiles of the Galactic halo and thick disk are listed in Table HI In Table HJ we also 
give the density profile of a hypothetically warped thick disk with the same warp properties as 
the thin disk. This profile is used in order to expl ore the possibility that the thick disk might 
be warped, as suggested by ICarraro et al.l (120081 ). Since the possible warp of the thick disk 
has not been explored in detail, the best assumption w e can make is that it is descri bed by 
the same parameters as the warp of the thin disk (e.g. lLopez-Corredoira et al.l 120021). This 



is a p lausible assumption since analytical models of warp formation (jSparke fc Casertano 
988 ) and N-body simulations, which model the formation of warped thin and thick disks 



Quinn et al.lll993l ) predict similar properties. 



The profile parameters, namely the RRLS density normalization p RR (for type ab 
RRLS), the power-law index n and flattening ratio c/a of the halo, as well as the scale 
length Iir and scale height h z of the thick disk, determined by several authors using differ- 
ent tracers are presented in Table \5 Since the only determination available to date of p RR 
for the thick disk is that of iLaydenl (119951 ). derived from nearby RRLS, we use this value 
(Pq R = 10 ± 4 RRLS/kpc 3 ) in the calculation of Nrr using the scale lengths and heights 
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derived by different authors (Table [5]). 



By solving Eq. [T] using the profiles given in Table HI we obtained the total number of 
RRLS expected in the range 3 kpc < D Q < 20 kpc. The results are shown in the last column 
of Table [5j in which Nr r has been com puted for RRLS of both types by correcting for the 
factor N a b +C /N a b = 1.29 (jLaydenlll995l ) to account for RRLS of type c. Also, we derived the 
number of RRLS as a function of heliocentric distance, which is shown in Fig. [6] together 
with the observed radial distribution of RRLS in our survey. From Fig. [SJ one can see 
that the theoretical RRLS radial distribution changes only slightly with the different sets of 
parameters. Most importantly, in the case of the warped thick disk there is a peak located 
at ~ 7.5 kpc, independent of the choice of the parameter set. This peak corresponds to the 
distance where the survey's line of sight intersects the highest density part of the thick disk's 
warp . The existen ce of this peak in the case of a warped disk is w hat led lMomany et al.l (120041 
20061 ) and others (lL6pez-Corredoirall2006l ; ICarraro et al.ll2008l ) to contest the interpretation 
of the CMa overdensity as a dwarf galaxy remnant, suggesting instead that it could be due 
to the warped Galactic disk(s). 

From Table [SJ one can see that 4 — 6 halo RRLS are expected in the distance range 
of the survey, and 2 — 3 or 3 — 6 RRLS are expected respectively from a normal or warped 
thick disk. This yields a total of 6 — 9 or 7 — 12 RRLS expected from a Galactic distribution. 
These numbers are in excellent agreement with the 10 RRLS found in our survey, leaving 
very little room (if any) for an excess of RRLS over the Galactic background. 

To take into account the random fluctuations that are obviously present in such small 
samples, we made 10 3 random realizations of a halo plus a normal thick disk (H+TD) and 
of a halo plus a warped thick disk (H+WTD). Each of these realizations was compared with 
the observed radial distribution, and the probabilities that the distributions are different 
were computed via a Kolmogorov-Smirnov test. These tests showed that for the H+TD 
and H+WTD distributions only 4% and 1%, respectively, of the 10 3 realizations were in- 
compatible with the observed distribution (i.e. P < 0.01). We conclude, therefore, that 
the observed RRLS radial distribution is consistent with the expectations from the Galactic 
models, although the small number of stars hinders any distinction between a normal and a 
warped thick disk. 

Several studies have shown that there is a correlation betwe en the |Fe/H| y alues of RRLS 
and their memberships in the stellar populations of the MW. iLaydenl (119951 ). for example, 
found that the majority of RRLS with [Fe/H]< —1.3 and [Fe/H]> —1.0 are members of the 
halo and thick disk populations, respectively, while the intermediate zone, —1.3 <[Fe/H]< 
—1.0, is occupied by a mixture of these populations. According to these limits, the probable 
memberships of the 8 RRLS in our sample that have [Fe/H] measurements are 4 halo stars, 
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1 thick disk star, and 3 in the mixture zone. While these assignments must be viewed with 
caution because the [Fe/H] distributions of both populations have tails that extend beyond 
the region of large overlap, they are consistent with the results obtained from the density 
profiles. 



4.2. Comparison with the Observed RRLS Radial Velocity Distribution 

As discussed above, we have measured the radial velocities of the 8 RRLS that lie in 
the distance range of the CMa overdensity (4 > D G > 12 kpc). In the upper diagram of 
FigJTl the radial velocity distributions of the RRab and RRc variables are depicted by the 
clear and the shaded histograms, respectively. They are plotted separately because the RRc 
sample is more prone to contamination from other types of stars. In the lower diagram, we 
have plotted, with arbitrary normalizations, the radial velocity distributi ons of A-F stars in 



the halo, thick, and thin disks, as given by the Besangon Galactic model (IRobin et al.ll2003l ) 
for the area of our survey. The halo and thick disk radial velocity distributions correspond 
to stars in the distance range of the RRLS in our survey; while the thin disk distribution was 
obtained for stars without any constraint on distance. The thin disk distribution is shown 
for reference since, as mentioned in §3.4[ some contamination may remain and these non- 
RRLS are most likely members of the thin disk. The dashed lines in the upper plot show the 
positions of the means of the distribution of each component. From this figure one can see 
that the observed radial velocity distribution is consistent with the expected halo and thick 
disk distributions. The radial velocities of some of the stars that we have classified as RRc 
are consistent with membership in the thin disk, but given the overlap in velocity with the 
halo and thick disk populations, it not certain that these stars are contaminants rather than 
true RRLS. In conclusion, it appears that 100% of the sample of RRLS can be explained 
by the contributions from the halo and thick populations and possible contaminants, and 
therefore, there is no excess of RRLS that requires another explanation. 



Martin et al.l (120051 ) measured the radial velocities of large samples of red clump (RC) 



and red giant branch (RGB) stars in the CMa overdens ity, and found that they define a 



sequence in a plot of D & against V r (see Figs. 6 & 12 in iMartin et al.l (120051 ) and our Fig. 
8). Their sample of RGB stars were observed in several two-degree fields that were centered 
on I = 240° and — 4?8 > b > — 10?8, and most of their sample of RC stars were observed in a 
field at (1, b) = (240?0, — 8?8). Our RRLS survey encompasses the field where the majority 
of RC stars were observed, and it overlaps with many of the fields where the RGB stars were 
observed (see Fig. 1). 



In Figure [HI the sequence that is defined by the observations of IMartin et al.l (120051 ) is 
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compared to our RRL S observations . For th e measurements of the RC stars, we adopted the 
mean values of V r that iMartin et al. (b()05h list in their Table 2 for different bins in D & . In 
our plot of these values in Figure El the vertical error bars depict the size of the bin in D Q 
(1 kpc for the RC stars), whil e the horizontal error bars depict the standard deviation of the 
Gaussians that IMartin et al.l fit to t he velocity pea ks. For the RGB stars, we adopted the 
mean V r and standard deviation that IMartin et al.l give for the distance bin of 7.5 — 11 kpc. 



For th e mean distance of these RGB stars, we adopted 8.5 kpc from Fig. 12 of IMartin et al. 
(bo05h . 



Figure El shows that several of the RRLS are widely scattered over the V r — D Q plane 
and that there is loose group at (V r , D Q ) ~ (70 km s _1 , 10.5 kpc). This group should not 
given much weight because 3 of the 4 stars are type c, the type most likely contaminated 
by thin disk interlopers, which would have instead D & J$ 4 kpc. Figure El s hows that the 
RRLS are not concentrated on the sequence defined by the observations of IMartin et al. 
J2005h . Two RRLS, the type ab variable 882 at (V r , D Q ) = (64 km s~\ 6.3 kpc) and the 
type c variable 39928 at 102 km s -1 , 10.7 kpc , lie w ithin the velocity and distance ranges 
of the RC and RGB stars that IMartin et al.l (120051 ) attribute to the CMa dwarf galaxy. 
While this could mean that these RRLS belong to this system, this is not necessarily the 
case. Members of the Galactic halo population are expected to have a very broad peak in 
velocity near 200 km s _1 and a roughly uniform distribution over the range of distances in 
Figure El (see Figuress 6 & 7). Consequently, every RRLS in Figure El is a candidate halo star. 
The behavior of the Galactic thick disk populat ion is depicted by the grayscale in Figure El 
which was obtained from the Besangon model (jRobin et al.ll2003l ). Since 882 and 39928 lie 
within the region expected to contain thick disk stars, it is possible that they belong to this 
Galactic population or the halo population and not to the CMa dwarf galaxy. Our survey 
has detected to 2 RRLS that could be members of the CMa dSph galaxy 



5. THE NUMBER OF RRLS IN DWARF SPHEROIDAL GALAXIES 



Out of the 10 RRLS found by our survey, at most 2 could be associated with the CMa 
dwarf galaxy. This is many fewer RRLS than expected. 

In Fig. [9j the number N RR of RRLS in dSph satellite gala xies of t he MW and M31 is 
plotte d vers us their absolute v isual magnitude My fdata from ISieeell (20061). iKuehn et al. 



( 2008 ) and Greco et al. ( 2008 ). and the compilation by Vivas &: Zinn (j2006 )). There 



is 



a strong correlation between the observed number of RRLS in a galaxy and its absolute 
magnitude in the sense that the number of RRLS grows with increasing luminosity. This 
correlation even holds for the extremely low luminosity satellites Boo, CVn I and II (My > 
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—8.5). It also holds for the Sgr dwarf which is at present undergoing tidal disruption. If 
the CMa dSph galaxy follows the same trend and has the esti mated My ~ —14.5 (estimates 
of the absolute V m agnitude range from —14.4 to —14.5, iMartmez-Delgado et al.l 120051 ; 
Bellazzini et al.ll2006l ). then it is expected to contain ~ 2880 RRLS. Our survey covers ~ 3% 
of the galaxy, and therefore, it should have detected ~ 86 RRLS. If CMa is more elongated 
with aspect ratio of ~ 6 : 1, the coverage of our survey would be > 2% and hence > 58 RRLS 
should have been detected by our survey. Since the expected number of RRLS is a factor of 
> 40 — 30 times greater than the number observed, we can safely reject the hypothesis that 
the CMa overdensity is caused by a "normal" dSph galaxy with My ~ —14.5. We can also 
firmly reject the possibility that the luminosity of this galaxy has been simply overestimated. 
Note that the trend in Fig. [9] indicates that ~ 600 RRLS are expected in a "normal" dSph 
with My ~ — 12, which is much fainter than any estimate for the CMa galaxy. Nonetheless, 
if such a system were present, our survey should have detected ~ 12 RRLS as an excess over 
the expected Galactic background of RRLS, which is not observed. 

It also clear from our survey that the putative CMa dSph galaxy must be unlike any 
other dSph galaxy in its production of RRLS per unit luminosity. Fig. [9] shows that the 
dSph galaxy that deviates by the largest amount below the line is the MW satellite Leo I. 
This was long considered to be an example of a galaxy whose stellar population is completely 
dominated by an inter mediate-age population because the re was no clear-cut evid ence for an 



older population (e.g. iLee et al.lll993l ; iGallart et al.lll999l ). The recent survey by iHeld et al. 



(120011 ) . which discovered 74 RRLS in Leo I, showed that even this exceptional system has 
an old stellar population containing RRLS. The Leo I galaxy has the lowest specific fre- 
quency^ of RRLS (S rr=1.9) of any of the 13 dSph galaxies that are brighter than My = —8 
(IVivas fc Zinnll2006l ). If a dSph galaxy of My = —14 produces the CMa overdensity, the 
results of our survey suggest that it has Srr < 0.2, an order of magnitude smaller than that 
of Leo I. 



6. DISCUSSION AND CONCLUSIONS 

The QUEST RR Lyrae survey in the CMa overdensity has yielded the following results: 



Ten RRLS, 6 of type ab and 4 of type c, were found by our survey, which spans ~ 20 
sq. deg. and is centered approximately at (/, b) = (240°, —9.5°). The distance range 



3 The specif frecuency, S rr,, is defined as the number of RRLS per unit absolute visual magnitude (My) 
normalized to My — —7.5 (jSuntzeff et al.lll99ll ). 
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covered by the RRLS is 3 kpc < D & < 20 kpc, and 8 out of the 10 RRLS are within 
the distance range of CMa overdensity (D@ < 11 kpc). 

• Integration of the density profiles of the Galactic halo and thick disk indicates that 
within the surveyed region the Galactic halo and thick disk should contribute a total 
of 6 — 9 RRLS, and if the thick disk is warped like the thin disk, the expected number 
is 7 — 12 RRLS. These numbers are in good agreement with the 10 RRLS found by 
our survey; consequently there is no excess of RRLS in the CMa overdensity. 

• The radial velocities of 8 out of the 10 RRLS were measured, and their values are 
consistent with membership in either the Galactic halo or the thick disk populations. 

• An analysis of the number of RRLS found in the dSph satellites of the MW and M31 
indicates that our survey should have discovered ~ 70 RRLS if the CMa overdensity 
is produced by a typical dSph galaxy with My = —14.5, the estimated luminosity of 
the CMa system. Even if this galaxy is 2.5 magnitudes fainter, which is much lower 
than any estimate in the literature, the expected excess of RRLS is ~ 12, which is 
not observed. The CMa dSph galaxy must have a significantly lower value of the 
specific frequency of RRLS (Srr) than Leo I, the Local Group dSph that has the 
lowest measured value. 



The hypothesis that the CMa overdensity is produced by a dwarf galaxy is severely 
constrained by the above results. We have shown that if it is a dSph galaxy, the type favored 
by the proponents of the hypothesis, then it must be unlike any of the satellite dSph galaxies 
of the MW and M31. No one has proposed that the putative galaxy is a dwarf irregular 



becau se there is no sign of a very young stellar population (age < 100 Myr, e.g. Ide Jong et al. 



20071 ) or star forming regions associated with the CMa overdensity. Nonetheless, we note 
that RRLS have been found in every dwarf irregular galaxy of the Local Group that has 
been adequately searched, an d that there is evidence for very old and metal-poor st ellar 
populations in the others (see IClementini et al.ll2003l ; iMomany et al.ll2005l ; lMateolll998l . and 
references therein). These populations have probably produced RRLS that await detection. 



According to Ide Jong et al.l (120071 ). the old stars in the CMD of the CMa overdensity 
have a mean metallicity of [Fe/H]~ —1 and ages in the range 3-6 Gyrs, although they could 
not constrain the presence of an older population. Clearly, the metals in this intermediate- 
age population must have come from earlier generations of more metal poor stars. There 
is no sign of a blue horizontal branch (BHB) in the CMD, and our survey has shown that 
there are very few, if any, horizontal branch stars in the instability strip. The remaining 
alternative is that the majority of the oldest and most metal poor HB stars of the CMa 
galaxy must lie on the red side of the instability strip and are therefore part of the red clump 
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that is observed in the CMD. However, at low metallicities (i.e. [Fe/H]< —2) a relatively 
young age is required t o keep the horiz ontal branch to red side of the instability strip (see, 
for example, Fig. 9 in iRey et al.ll200ll ). Thus, the CMa dwarf must have started its star 
formation at a later time than is observed in other dwarf galaxies. This was once thought 
to be the case for Leo I, but we now know that this was incorrect ( jHeld et al.ll200ll ). 



Martin et all fl2004af ) have suggested that the globular clusters NGC 1851, 1904, 2298, 



and 2808 once belonged to the CMa dwarf galaxy. These clusters span a range in horizontal 
branch type from predominantly blue to p redominantly red, although each o ne contains 
many BHB stars and at least a few RRLS (ILee et al.lll994l ; iForbes et al.ll2004l ). The idea 
that they originated in the CMa dwarf is hard to reconcile with the lack of a BHB in the 
CMa CMD and our result that there is no excess of RRLS in the CMa overdensity. In other 
dwarf galaxies that have globular cluster systems (e.g., the Sgr and Fornax dSph galaxies), 
there are stellar populations in the main body that resemble the ones in the retinues of 
globular clusters. 

In summary, the dwarf galaxy explanation for the CMa overdensity requires that the 
system be unique among dSph galaxies, and perhaps also dwarf irregulars, in its low pro- 
duction of RRLS. The absence of RRLS in the overdensity is consistent with the alternative 
viewpoint that a combination of the t hin and thick disk and sp i ral arm populations of 



the MW produces the CMa oyerdensitv (jMomany et al 



Lopez- Corredoira et all 120071 ; ICarraro et al. 



2005 



2008 



2004 



2006 



Lopez- Corredoira 



Moitinho et al 



2006 



2006; 



Powell et al. 



20081 ). However, it does not prove that this alternative is correct, and more research on this 
question is warranted. 
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Fig. 1. — Schematic plot of the areas in CMa studied by different authors. The blue area 
shows the QUEST RRLS survey footprint. The colored rectangles represent the areas studied 
by Bellazzini et al 2006 (purple), Martinez-Delgado et al. 2005 (red), Carraro et al. 2005 
and Moitinho et al. 2006 (cyan), Butler et al. 2007 and de Jong et al. 2007 (yellow), Conn 
et al. 2007 (magenta), Powell et al. 2008 (dark pink) and Carraro et al. 2008 (green). The 
axes of the dashed-lined ellipse are equal to the full-width at half maxima derived for CMa 
usi ng an aspect r atio o f 5 : 1 and the fu ll-width at ha l f max imum in the b direction measured 
by iButler et al.l (120071 ) . As noted by iButler et al.l (120071 ). the aspect ratio of CMa could 
be larger than 5 : 1, in which case the ellipse would be more elongat ed in longitud e . The 
colorscale represents the Ay values obtained from interpolation of the ISchlegel et al.l (119981 ) 
dust maps. 
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Fig. 2. — Typical standard deviation of PSF magnitudes versus mean magnitude, based in 
~ 20 observations of ~ 80000 stars. Solid triangles represent the positions of the RRLS 
identified in this work. 
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Fig. 3. — The light curves of the identified RRLS. Observations made with the QUEST 
camera are shown as open triangles, and follow-up observations from the SMARTS telescopes 
and the NOV 1.0m Reflector are shown with open squares and circles respectively. Solid 
lines show the best fitting template. Typical error bars are smaller than the symbol size. 
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Fig. 4. — The light curves of the identified W UMa eclipsing contact binaries. Symbols are 
as in Fig. [31 Solid lines show the best fitting template. 
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Fig. 5. — Completeness of the survey as a function of V magnitude. 
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Fig. 6. — Observed and expected numbers of RRLS as a function of heliocentric distance. 
From top to bottom: Observed distribution, Galactic halo, normal thick disk and warped 
thick disk distributions. 
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Fig. 7. — Upper panel: Radial velocity distribution of the 8 RRLS with V < 16.5. The shaded 
histogram corresponds to the radial velocity distribution of RRc stars. Bottom panel: Radial 
velocity distributions of thin disk A-F stars (dashed dotted line), thick disk (dashed line) 
and halo (solid lin e) stars in the surveyed area, taken from the Besangon Galactic model 
( iRobin et al.l 120031 1 . The dashed lines in the upper plot show the positions of the means of 
the distribution of each Galactic component. 
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Fig. 8. — Heliocentric distance versus radial velocity of the 8 RRLS with V < 16.5. Solid 
triangles represent RRLS of type ab (up) and c (down). Open symbols represent RC (rhom- 
buses) and RGB (pentagon) stars identified with the CMa overdensity Martin et al.l (120051 ) . 
The grayscale shows th e distribution of th ick disk stars in the surveyed area from the Be- 
sangon Galactic model (IRobin et al.ll2003l ). The distribution of halo stars in this diagram 
(not shown) has a broad peak in Vr at ~ 200 (see Fig. 8) and gentle falloff with D Q . 
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Fig. 9. — The observed number of RRLS in dSph satellites of the MW (circles) and M31 
(triangles). A least squares fit to the data is represented by the strai ght line. Note that 
> 1000 RRLS are expected in the CMa dSph galaxy if it has M v 14 jMartin et alil2004al : 



Bellazzini et al.l 120041 ; iMartmez-Delgado et al.l 120051 ) and follows the trend in this diagram. 
The position of Leo I, the MW dSph galaxy with the lowest Srr (see §4), is marked. 



Table 1. Summary of previous stellar population studies conducted in CMa. 



Tracer 


Age (Gyrs) 


[Fe/H] 


D& (kpc) 


Technique 


Filters 


Interpretation 


Reference 



Blue Plume 
Blue Plume 
Blue Plume 
Blue Plume 
B5-A0 stars 


1-2 

- 0.7 - 2 

<o.i 
<o.i 


-0.3 or 1.0 

-0.37 to -0.5 


~ 7.5 or 9.3 
10.8 

6.0 ±2.7 

9 8+ 1 ' 5 
a,o -1.0 


MS-fitting 
HD modelling 
TCDs, MS-fitting 
TCDs, spectroscopy 
TCDs 


BV 
BR 

UBVRI 

UBV 

UBVRI 


dSph galaxy 
dSph galaxy 
Spiral arm 
Spiral arm 
Spiral arm 


B04 
dJ07 

C05, Moi06 

P08 

C08 


MS 


3-6 


~ -1.0 ±0.1 


~ 7.5 


HD modelling 


BR 


dSph galaxy 


dJ07, B07 


RGB, MS 


4-10 


-0.66 to -0.35 


~ 8.1 ± 1.2 


Isocrone fitting 


BV 


dSph galaxy 


B04 


M-giants 


4-10 


-1.0 < [M/H] 


7.1 ±0.1 


CMD fitting 


JHK 


dSph galaxy 


M04 


RGB, RC 






5.5 to 8.5 


Spectroscopy 




dSph galaxy 


M05 


F-G-K stars 


6±2 


0.3 ±0.3 


6 


TCDs 


UBVRI 


Warped Thick Disk 


C08 



CO 



Refer ences. — B04 = Bellazzini et alj ll2004h. M04=lMartin et alj J2004ah. M05=lMartin et alj lipoih. C05=ICarraro et al.l ll2005l), 
Moi06 = lMoitinho et all (120061 '). B07= lButler et alj ll2007t ), dJ07= lde Jong et alj faOOTj) , C08= ICarraro et alj fcOOSf) . P08= |Powell et alj 

(2008) 



Table 2. Light curve and physical parameters of RRLS. 



ID 


a 


8 


N 


Period 


Ampv 


Type 


H J Dmax 


(Vb> 


D Q 




aV~f 


[Fe/H] 




(deg) 


(deg) 




(d) 






2450000 d 




(kpc) 


(km s- 1 ) 


(km s- 1 ) 




7701 


108.429950 


-28.128540 


16 


0.58446 


0.71 


RRaft 


3260.75930 


14.22 


5.4 


200 


15 


-1.51 


882 


104.492882 


-27.708120 


29 


0.46353 


1.18 


RRaft 


3372.89115 


14.54 


6.3 


64 


15 


-1.52 


20594 


105.830994 


-27.565201 


35 


0.50893 


1.32 


RRafe 


3372.41209 


15.73 


10.9 


72 


16 


-1.83 


26896 


108.468697 


-29.726669 


29 


0.39725 


0.96 


RRafe 


3372.85468 


13.57 


4.0 


8 


15 


-1.29 


29821 


108.591827 


-31.073851 


25 


0.37719 


0.25 


RRc 


3372.83819 


15.33 


9.0 


54 


15 


-1.27 


31038 


107.737091 


-27.722931 


25 


0.26672 


0.31 


RRc 


3372.95514 


16.73 


17.2 








32697 


107.872360 


-27.761780 


25 


0.20140 


0.46 


RRc 


3372.81509 


15.54 


10.0 


61 


14 


-1.26 


34757 


109.645363 


-29.783131 


20 


0.54532 


1.20 


RRafe 


3372.59181 


17.09 


20.3 








39928 


109.702118 


-30.578770 


22 


0.33740 


0.33 


RRc 


3372.86688 


15.69 


10.7 


102 


21 


-0.78 


40676 


107.920067 


-28.308350 


31 


0.57517 


0.80 


RRab 


3373.13496 


15.54 


9.9 


182 


15 


-1.44 



Table 3. Light curve parameters of rejected RRLS candidates. 



ID 


a 


5 


N 


Period 


Vamp 


H J Drnax 


(Vo) 


Reason for 




(deg) 


(deg) 




(d) 




2450000 d 




Rejection 


5564 


105.416870 


-28.277680 


22 


0.29700 


0.17 


3258.68871 


15.69 


latc-F spectra 


10651 


105.439903 


-31.074100 


30 


0.37339 


0.37 


3372.46691 


13.88 


W UMa 


25711 


105.681969 


-28.129801 


32 


0.42189 


0.42 


3372.97294 


15.03 


W UMa 


6460 


107.042080 


-27.251020 


22 


0.21343 


0.29 


3259.82428 


14.49 


late-F spectra 


18056 


107.396880 


-28.883540 


19 


0.27842 


0.17 


3259.93843 


14.10 


W UMa 


17237 


108.275560 


-26.897120 


17 


0.24430 


0.24 


3259.75716 


14.36 


late-F spectra 
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Table 4. Density profiles of the Galactic halo and thick disk. 



p(f) 



Halo p RR^ X 2 + y 2 + { _*_ ) 2 



Thick Disk p* R exp(- ga ^ K0 ) exp(-g) 



Warped p* R exp(- ^~ fl0 ) exp(-^^) 

Thick Disk with z w =[c w R gal (pcy™ s in(^-^)-i5]pc if i? 9a! <i3k P c 

With z w =z w (13kpc) if fi 9a! >13kpc 



Note. - In the table, x, y, z are cartesian Galactocentric 
coordinates, R & = 8 kpc is the distance from the Sun to the 
Galactic center and R ga i the Galactocentric distance along the 
plane. In the density profile assumed for a warped thick disk, 



5.25 and 



-5°, are th e warp 



C w = 2.1 x 10" 19 pc, e v 
parameters according to lLopez-Corredoira et al.l (120021 ) . 





Table 5. 


Numb 


er of RRLS expected 


in each Galactic component. 






Reference 




RR (#.) 


n 


(c/o) 


Tracer 


Nrr 


Halo 


Vivas & Zinn (2006) 
Wetterer & McGraw 1 
Preston et al. (1991) 


19961 


4 O+0.5 
4 - z -0.4 
5.3 ±2 
4.8 


-3.1 ±0.1 
-3.53 ±0.08 
-3.2 ±0.1 


Variable 
Variable 
Variable 


RRLS 
RRLS 
RRLS 


6 
4 
6 




Suntzcff ct al. (1991) 




4.5 ±1 


-3.5 


1 


RRLS 


5 




Reference 




RR I # 1 

Pq ( J^> 


h z (pc) 


h R (kpc) 


Tracer ff™ 1 N™ R rped 


Thick 
Disk 


Brown et al. (2008) 
Cabrera-Lavers et al. 
Larscn & Humphreys 
Oiha (2001) 


(2005) 
(2003) 1 " 




1026 ± 100 
1061 ± 52 
900 ± 65 
860 ± 200 


3.5 ± 0.5 
3.04 ±0.11 
4.7 ±0.2 

q y+0.8 
d ' ' -0.5 


BHB stars 3 
2MASS RC stars 2 
Star counts 3 
2MASS star counts 2 


5 
3 
6 
4 



Note. — The symbol Nrr represents the number of RRLS in the distance range from 3 to 20 kpc, both for the halo and normal/warped 

thick disk. '''For the scale-height Larsen & Humphreys (2003) obtained hz = 870 — 930 ± 50 — 80pc, for the present 
study we use the intermediate value reported here. 



